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' ^ , We analyse the measured critical current I m in a mesoscopic 4-terminal S/N/S structure. The 

pr^ ■ current through the S/N interface is shown to consist not only of the Josephson component I c sirup, 

but also a phase-coherent part I sg cos ip of the subgap current. The current 7 m is determined by 
the both components I c and I sg , and depends in a nonmonotonic way on the voltage V between 
superconductors and normal reservoirs reaching a maximum at V = A/e. The obtained theoretical 
resultas are in qualitative agreement with recent experimental data. 
Recent achievements in nanotechnology have revived interest in the study of nonequilibrium and phase-coherent 

q . phenomena in superconductor- normal metal (S/N) structures. One of the most remarkable, discovered recently 
0, was the observation of the sign reversal of the Josephson critical current I c (the so-called 7r-junction) in a multi- 
terminal mesoscopic Nb/Au/Nb structure under nonequilibrium conditions. By passing an additional current through 
the N layer or, in another words, by applying a voltage V to the normal reservoirs (see Fig.l) with respect to the 
superconductors, one can create a nonequlibrium electron-hole distribution, or at least one can shift this distribution 
with respect to the electron-hole distribution in the superconductors. Under this condition, the critical current I c 
decreases with V and changes sign at a certain value of the applied voltage V. This effect was predicted first in Ref. 
where a ballistic 3-terminal structure was considered (for more details, see also Refs. [0H). In diffusive 4-terminal 
S/N/S structures, the sign-reversal effect has been considered in Refs. ®^7$ (see also The sign-reversal effect 

and switching of the 7r-junction into a state where ip = n has much in common with an instability of an uniform 
superconductor with a nonequilibrium distribution function |lO| , |li"| ] . 

In multi-terminal S/N/S structures one can observe not only the sign reversal effect, but also a number of other 
interesting phenomena. For example, the conductance of a normal wire between N reservoirs oscillates with varying 
phase difference tp (see review articles 12 In addition, as shown in Refs. P,|l4||, the measured critical current /„ 
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depends on the geometry of a particular structure and instead of decreasing may also increase with increasing voltage 
V. In particular one can observe Josephson-like effects (plateau on the h(Vs) curve, oscillations of the measured 
C\| . critical current I m in a magnetic field etc) even if the Josephson coupling between superconductors under equilibrium 
conditions is negligable. The reason for these effects is that the current I m in a multi-terminal S/N/S structure is 
determined not only by the Josephson component I c shx(p, but also by the phase-dependent subgap current I sg cos ip 
through the S/N interface. Therefore even in the case of a small I c , the current I m can be altered by varying the 
phase (j>. An increase of the critical current was observed in the recent paper [fl5| where a mesoscopic three-terminal 
S/N/S structure was studied. The authors used a third superconductor as a reservoir the electric potential of which 
was shifted with respect to the other two superconductors by the voltage V. The measured critical current reached its 
maximal value when the magnitude of V was comparable with A. At some, not too low temperatures T the measured 
critical current I m exceeds its magnitude in the equilibrium state: I m (V) > J m (0). In the present paper we show that 
the enhancement of the supercurrent observed in Ref. |15| is most likely caused by the mechanism mentioned above. 
q , In Refs. ||,|l4| the model case of gapless superconductors was considered where there is no singularity in the density-of 
states in superconductors at e = A. Here we will consider the case of ordinary superconductors with an energy gap A 
and show that the enhancement of the critical current reaches a maximum for V of order A. The voltage dependence 



Im (V) calculated for different temperatures is in qualitative agreement with the experimental data. 



We consider the structure shown in Fig.l which differs from the structure studied experimentally. However in 
our opinion, this difference is not essential and allows us to give at least qualitative explanation for the phenomena 
observed in Ref. Jl5[ . First, we assume for simplicity that the structure under consideration is symmetrical, i.e. it 
has four terminals and not three as in the experiment. Secondly, we consider normal reservoirs in order to avoid 
complications which would arise in case of superconducting reservoirs (ac Josephson effects when the finite voltage is 
applied to the S reservoir). We also assume for simplicity that the contacts between the N wire and N reservoirs are 
good (the resistance of the N wire/N reservoir interface is much smaller than the resistance of the N wire), whereas 
the S/N interface resistance is finite (larger or less than the resistance of the N wire). We will study the diffusive case 
which corresponds to the experiment fll5|]. 
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In order to find the dependence of the effective critical current I m (V) ( the definition of I m (V) will be given 
later), we need to determine two distribution function / + and /_. Both these functions are isotropic in space. 
The function / + is related to a symmetrical part of the distribution function in the electron-hole space:/+(e) = 
1 — (n-f(e) +pi(e)) — 1 — (ftj_(e) + Pt( £ ))> here Pl( e ) = 1 — n l(~ £ ) is the hole distribution function. It determines the 
critical current I c . The function / _ describes the electron- hole imbalance and determines the electric potential and 
current: /(e) = — (n-f(e) — Px( e )) = — ( n l( e ) ~ Pt( e ))- Equations for / + and /_ are obtained from an equation for the 

matrix Keldysh function G (see, for example For the structure shown in Fig.l they can be written in the form 

Ld x (M_d x f-(x)+J s f + -J an d x f+(x))=r[A_5(x-L 1 )+A- 5(x + Li)]. (1) 

Ld x {M + dJ + {x)+J s f_+J an dJ_{x)) = r[A + 6(x - L,)+ A+ S(x + I*)] (2) 

A A 

Here all the coefficients are expressed through the retarded (advanced) Green's functions G R — G R a z + F R and are 

A A 

equal to M± = (1 - G R G A T [F R F A ) x )/2; 

J an =(F R F A ) z /2, J s = (l/2)(F R d x F R - F A d x F A ) z , A_ = {vv s + g 1+ )f_ - (g z -f eq + g z+ f + ); A + = {vv + 

<?!-)(/+ f eq ) - 9z-f-;9i± = (1/4)[(F R ±F A )(F R ± F A )] X ; g z± = (l/4)[(F R T F A )(F R ± F A )] Z ; 

The coefficient r = R/Rb, R = pL/d is the resistance of the N film per unit length in the z-direction,p is the specific 
resistivity of the N film, d is the thickness of the N film, Rb is the S/N interface resistance; the functions A- and A+ 

A A 

coincide with A_ and A + if we make a substitution <p — > — <p. We introduced above the following notations (F R F A )i = 

A A A A A A 

Tr(F R F A )/2, (F R F A ) Z = Tr(a z F R F A )/2 etc.; v and v s are the density-of states in the N film at x = L x and in the 
superconductors. The boundary conditions for / + and /_ are: /+ (L) = Fy+ and f-(L) = Fy-', the functions Fy± 
are the corresponding distribution functions in the normal reservoirs: Fy± — [tanh((e + eV)(3) ± tanh((e — eV)j3)]/2. 
We set the electrical potential at the superconductors equal to zero and assumed that the width of the S /N interfaces 
w is small compared to L\^- 

Eq.(l) describes the conservation of the electric current (at a given energy). The term in the brackets on the left is 
the total partial current in the N wire, consisting of the quasiparticle current (the first term), the supercurrent in the 
interval (— L\,Li) (the second term) and a " nonequilibrium supercurrent" (the third term). The coefficient M is a 
quantity which is proportional to the diffusion coefficient renormalised due to proximity effect. The right hand side is 
the partial current through the S/N interface; the term (vvs + <7i+)/- is the quasiparticle current above (vvsf-) and 
below (<7i + /_) the gap. The term (g z _f eq + g z +f+) is the Josephson current in nonequilibrium conditions. Eq.(2) 
describes the conservation of the energy flux (at a given energy). The coefficient A + is zero below the gap ( complete 
Andreev reflection) as the difference (F R — F A ) equals zero at e < A. 

The solutions of Eqs.(l)-(2) can be found exactly and expressed in terms of the retarded (advanced) Green's 
functions which obey the Usadel equation. First we note that the expressions in brackets in the left hand side of 
Eqs.(l)-(2) in the regions (0,Li) and (L\,L) are equal to the constants of integration Ci,2±- The constants C\ t 2- 
relate to partial currents (Ci,2- = eJi^p/d). The partial currents Ji ; 2 are the currents per unit energy and 
connected with the electrical currents via the relation 

/>oo 

h.2 = / deJi, 2 (e) (3) 
Jo 

Our aim is to find the current J3 and express it in terms of the control current I2 (or voltage V) and the phase 
difference ip. We note that the distribution functions f±(x) are constants in the region x £ (0, L\) and vary in the 
region x G [L\, L) reaching Fy± at x = L. Dropping details of calculations, we present final results for limiting cases. 

a) Large interface resistance: r « 1. 

One can show that in this case /+(0) = {Fy+ + feq{r2vv s ))/(l + r 2 vv s ) and /_(0) = Fy-/(1 + T2W S ), where 
T2 = r(L-2/L). The current ^3 through the S/N interface consists of three terms 

h(V) = I {V) -I C (V) sin <p + I sg (V) cos (4) 

Two of them (/ , I sg cos ip) are the quasiparticle currents and one (I c sin ip) is the Josephson current. This expression 
shows that at a given control voltage V and zero voltage difference between the superconductors [if is constant in 



2 



time) the current I3 may vary with changing ip in the limits:) h{V) — I (V) |< I m (V). This means a plateau on the 
Vs(h) characteristics (see HQ); here Vs — (h/2e)d t p is the voltage difference between superconductors. We can 
write the phase-dependent part of I3 in the form = I m sin(<p + a), where J m = \Jl% + Ig g is the measured critical 

current, cos a = —I c /I m . In the considered limit of high interface resistance, we have for I c and I sg 
I c {eR b ) = - / °° de{lm(F s {F y - F x ))f o (0) + ReF s Im(F y - F x )/(0)} 
I sg (eR b ) = / °° deg sg f(0) = / °° de lmF s lm(F y - F x )f(0); 

Here 8 = k e L,02 = k e L>2,k e = <J (~2ie + j)/D, 7 and D is the damping rate and diffusion coefficient in the N 
film, g sg = 9i + is the normalised subgap conductance (see the expression for A). The functions F yi F x are the 



components of the retarded Green's function in the N film: F R = F x io x + F y ia y , and F$ = A/ y/{e + iT) 2 — A 2 
is the amplitude of the retarded Green's function in the superconductors. If we linearise the Usadel equation, we 
obtain F y — F x = 2F$ sinh 2 62 / (0 sinh 29) . We note that the numerical solution of the Usadel equation shows that the 
linearised solution is a good approximation even if r = 1 (at r = 1 the difference between the exact and the linearized 
solutions at the characteristic energy e = el = D/L 2 is less than 5 percent). In Fig. 2 we plot the V dependence of 
I c , I sg and I m where we see that the real critical current I c decreases and changes sign with increasing V , whereas the 
measured critical current I m first decreases and then increases. Its maximum may exceed I c (0). The reason for such 
a behaviour of I m is the third term on the right side in Eq. (5) which describes a contribution of the phase-dependent 
part of the subgap quasiparticle current I sg through the S/N interface to the current I3. The current I sg is zero at 
V = and increases with V; this current leads to a low |l6| and high temperature peak in the conductance. Its 



phase dependence was measured in Ref. |18 and discussed in many papers (see review articles |l2j,|l3j). One can see 
from Fig. 2 that due to the current I ag the measured critical current I rn remains finite when I C {V) turns to zero. 

Fig. 3 shows the dependence of the measured critical current I m on the control voltage V for different temperatures. 
Our results qualitatively agree with the experimental data of Ref. ]l5| | ; that is, the current I m reaches a maximum at 
V = A/e and this maximum exceeds the equilibrium value of I c — I c (0) when the temperature is not too low. One 
can see, in agreement with the experimental results of Ref. ]I5| ] , the maximal value of I m depends on the temperature 
much weaker than I c . Although it is difficult to carry out a quantitative comparision between theory and experiment 
because in the experiment the width w and the interface resistance Rb were comparable with and R respectively, 
and a superconducting reservoir was used instead of a normal one (therefore, strictly speaking, one must take into 
account ac Josephson effects) . 

An important point to note is that our results do not mean that the sign reversal of the real critical current I c can 
not be identified directly. Consider for example a fork-shape circuit; this means that two vertical superconducting 
leads in Fig.l are attached to a T-shape (inverted) superconducting lead. Analysing the stability of the state with 
negative I Cl one can easily show that the state with ip = is unstable with respect to fluctuations of tp and the 
system switches to a state with a circulating current. Indeed, taking into account the fluctuating voltage at the 
superconductor Vs = hdtp/2e, we replace V in Eq.(4) by V — Vs- We then write down the equation for the current ^3 
in the lead attached to the left superconductor; this equation coincides with Eq.(4) if p> is replaced by —p. Subtracting 
these equations for ^3 and J 3 , we arrive at the equation for a circulating current 7 c j r = — (I3— ia)/2: 

Idr = I C {V) sin ip + Vs(Ro + Rsgcosp). (5) 

where i?o = dI /dV and R sg = dI sg /dV . Fluctuations of J c j r lead to a magnetic flux $ = I C i r L/c in the loop 
which is related to p: $ = ^o'P, here $ G is the magnetic flux quantum and we assumed the absence of flux in the 
ground state. We find readily from Eq.(5) that the state with p = is unstable if I C (V) < and | I C (V) \> c& /L, 
where L is the loop inductance 

b) Small interface resistance. 

One can show that in this case the function /_ (0) is zero in the main approximation with respect to the parameter 
(r9)~ 1 (this means that the condition r 2 >> A/cl should be satisfied; here cl = D/L 2 is the Thouless energy). The 
function / + , which determines the Josephson current, in the main approximation is equal to Fy+ at e |< A and to 
f eq at I e |> A. Therefore the dependence I C (V) is similar to that found numerically in Ref. for another geometry 
(for small interface resistance); that is, the critical current I C (V) changes sign with increasing V at V of the order 
of the Thouless energy. As to the current I2, it does not depend on the phase difference in the main approximation. 
Indeed, in order to find I2 we need to solve the Usadel equation in the region x 6 (Li,L) with boundary condition 

A A A A A A 

which is reduced to G R =Gg. Making the gauge transformation G^=> S G R S + , we can exclude the phase (here 

A 

S= cos((^/2) + icr z sin((yj/2)). Therefore in the main approximation the third term in Eq.(4) is zero. 
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In conclusion, we have studied the dependence of the measured critical current I m on the voltage V between normal 
reservoirs and superconductors in a 4-terminal S/N mesoscopic structure. The current I m is shown to decrease with 
increasing V, then to increase reaching a maximum at V = A/e. Our results qualitatively agree with experimental 
data obtained in the recent paper fl5| . 

We are grateful to the EPSRC for their financial support. 
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FIG. 1. Schematic view of the 4-terminal S/N/S structure under consideration. The electric potential of the superconductors 
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FIG. 2. The measured (Im) and real (I c ) critical currents vs the control voltage V. The amplitude of the phase-dependent 
part (I sg ) of the subgap current is shown by the dashed line. The currents and voltage are measured in units CLR/eR 2 . and 
£L/e respectively ( e_L = hD/L 2 is the Thouless energy). The parameters are: A — 4e_L,T = ei/4,Li/i = 0.3, r — 0.3. 




FIG. 3. The measured critical current (Im) vs V for different temperatures: /3 = €l/2T. The parameters are: 
A = 10e L ,Li/L = 0.3, r = 0.3. 
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